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Abstract 

An Lxoo system of odd number of coupled Heisenberg spin chains is studied using a 
degenerate perturbation theory, where L is the number of coupled chains. An effective 
chain Hamiltonian is derived explicitly in terms of two spin half degrees of freedom 
of a closed chain of L sites, valid in the regime the inter-chain coupling is stronger 
than the intra-chain coupling. The spin gap has been calculated numerically using the 
effective Hamiltonian for L = 3, 5, 7, 9 for a finite chain up to ten sites. It is suggested 
that the ground state of the effective Hamiltonian is correlated, by examining various 
variational trial states for the effective-spin chain Hamiltonian. 
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There is a growing interest in coupled chains of Hubbard-Heisenberg spin systems, 
following the experimental realization of coupled arrays of metal-oxide-metal ladder systems 
3- A number of investigations f|, ^, |6|, |7| of weakly-coupled Heisenberg spin chains 
have been carried out, in which the coupling between the chains is weaker than intra-chain 
coupling, and have provided a strong indication that a system of even number of chains can 
be understood in terms of a short-range resonating valence bond (RVB) picture, with spin 
gap and a finite spin-spin correlation length. Whereas a system of odd number of chains 
is gapless with power-law spin correlations, indicating that it is in the same universality 
class as a single spin chain. The framework for understanding the difference between even 
and odd numbered chains is the RVB picture ||, |8|], which suggests that an even- numbered 
coupled chain system can be thought of as an integer spin chain exhibiting a Haldane spin 
gap, and an odd-numbered coupled system maps onto a half-odd integer spin chain. Below 
we will construct an effective spin Hamiltonian for odd-numbered coupled system, which is 
cast explcitily in terms of spin half degrees of freedom. This can be done straightforwardly 
in the regime the coupling along the chain is weaker than the coupling between the chains, 
complementary to the earlier studies. 

We consider the Heisenberg Hamiltonian given as 



defined on an L x N lattice where I = 1,2, .L and n = 1,2,..N, and s^ n is a spin-1/2 
operator at site (l,n). The intra-chain exchange coupling constant is J' and the interchain 
coupling constant is J, and h is a magnetic field. We use periodic boundary conditions 
in both directions, and as we see below it makes a big difference. We are interested in 
regime where J' < J, and the case of odd-numbered coupled chains, L = 3,5.. We can 
think of the system as consisting of N sites with 2 L degrees of freedom at each site, and a 
complicated bond interaction between the sites mediated by the J' bonds in the Hamiltonian 
given above.. Starting with an unperturbed system with J' = 0, we proceed to include the 
effect of a non-zero J' within a perturbation theory, This is achieved by identifying the most 
relevant degrees of freedom for a given site, and projecting out the not so important degrees 
of freedom. Below we explicit show an analytical construction of the effective-spin chain 
Hamiltonian for the case of L = 3, and 5, and for a larger L one can a similar procedure can 
be implemented numerically. 

Let us consider the case of J' = 0, and h = 0. The ground state of the system is 
direct product of N copies of the ground state of a closed chain with L sites. The ground 
state of Heisenberg Hamiltonian of a closed chain with an odd number of sites is in sector 
with total spin S = | and S z = ±|, implying at least a two-fold Kramers degeneracy. As 
a matter of fact the ground state is four-fold degenerate. This extra degeneracy arises from 
chirality. This can be viewed as an unpaired spin in an otherwise singlet ground state, i.e. 
a spinon which is forced to be there in the ground state 0. In addition to the spin, the 
spinon carries a finite ground state momentum ±4ir/L (orbital angular momentum about 
the center of the closed chain). In the case of an open chain the chiral degeneracy disappears, 
and we will be left with a two-fold Kramers-degenerate ground state with S z = ±1/2. The 





2 



two ground states in S z = \ sector can be written as \ , = a ±({ n i})\{ n i} >> where Hi 

2 ,± 

refer to locations of down spins {n\ < n 2 ...). Similarly the ground states in S z = — ~ sector 
cab be written (here refer to the up spins). The amplitudes a± are found using Bethe 
anstaz[[L0|; for L = 3 the amplitudes are given by (there is only one down spin in this case) 
a±(ni) = exp(±i2imi/3), and for L = 5 the amplitudes are given by (there are two down 
spins in this case) 

a ± (ni,n 2 ) = exp ±2ik(ni + n 2 ) cos((k - r]/5)(n 2 - ni) + T)/2), (2) 

where k = 2n/5 and the phase shift rj = 2.243... These two different states have spin 
currents going in two different directions. We refer to these two states as having chirality 
by defining x z \ ( t ) ± > = ^ll^i >J anc ^ ^ e assoc i a ted raising and lowering operators through 
X ± \± >— 0, x ± l = F >= |=t >. The chirality operator \ can be written in terms of the 
original spin operators. Following [pT| , we define a permutation operator Pi at a site / which 
permutes the spin labels such that / + m — > N + I — m. This is equivalent to doing a 
reflection on a regular polygon with L sites, around a straightline bisecting the angle at site 
/. The P operators can be readily constructed from the spin operators, as each of them 

involve (L — l)/2 pair-wise permutations, as Pi = 2 2 rim=i( s z+m-SAr+;_ m + ¥)■ Also, Pi 
can be expressed in terms of chiral lowering and raising operators as Pi = A^x + + \ 
where Ai = exp(i2nl / L) . The chirality operator is now related to the original spin operators 
through [Pi,Pi +m ] = 4ix z sin27rm/L, and x x = -Pi/2. 

If we denote the gap in the spectrum for this L— site closed chain by 8l (see Table 
1 for the actual values of the gap for L = 3,5, 7, 9), for J' < 5l it is a good approximation 
to just keep the ground state manifold, 4 Ar -fold degenerate for J 1 = 0, and dropping the 
excited states. This is expected to suffice for understanding the low-energy behaviour of the 
system described by the Hamiltonian given in Eq. ??, as the excitation processes to include 
the the states dropped have high energy. The ground state subspace we keep, belonging 
to spin-1/2 sector of each site (of 2 L states), is identified as the most relevant degrees of 
freedom. For L = 3, this procedure keeps all the spin-1/2 at each site. However, for larger 
L > 3, there are also spin-1/2 excited states with large energy which are eliminated in our 
procedure. This approximation is not so good for L = 9 or more, as the the range of J' for 
validity of our procedure becomes very restricted, viz. J' < 0.8612 for L = 9. The ground- 
state subspace for J' = 0, is only 2 N dimensional if we use open boundary conditions for 
the L— site chain, in contrast for periodic boundary conditions we get a much larger Hilbert 
space, rendering the perturbation theory to be valid in a smaller regime of J' as the gap 
for the L— site open chain is smaller than that of a closed chain. Each of the basis states 
-0j for the full system is a direct product of N ground state eigenfunctions of close chain of 
L sites, ipi = 0^ x <f)j 2 ... x <j)f , where Zj = 1,4. To find the effective Hamiltonian within 
this subspace, we compute the matrix elements < ipilH'lipj >, and use them to construct an 
effective Hamiltonian. The prime denotes that only J'bonds are included. The evaluation 
of the matrix elements reduces to a two-site problem (actually two closed chains of length L 
connected by L bonds of strength J') problem as we have pair-wise interactions only. 

Let us use two spin-^ effective degrees of freedom for constructing an effective Hamil- 
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tonian, the spin and the chiral quantum numbers Si and the chirality %jof the L— site closed 
chains. An explicit calculation shows that the effective interaction strength Aij between two 
closed chains connected by a J' bond (the corresponding operator is s^i-s*,-^) at sites i and j 
is given in terms of the total spins of the two closed chains and the permutation operators, 
which we discussed above, at sites i and j respectively of two chains as 

= J'a L S..S, ( lL + ^) ( 7£ + 5^)). (3) 

Here for L = 3, a 3 = 16/9, 73 = 1/4, and for L = 5 a 5 = 16/3, 75 = l/20y/p. The parameter 
P is given in terms of x 1 = cos k — rj/5 + rj/2, and x 2 = cos2(A; — rj/5) + rj/2 as 

2 (1 + cos 2fc)xf + (1 + cos k)x\ — x\x\ , . 

5 x\ + x\ 

For L > 5, the wave-function renormalization constants and 7i are difficult to calculate 
analytically, as the number of s 2 -basis states in S z = 1/2 becomes very large. However, 
the constants can be calculated numerically, and the form of the interaction remains valid. 
Numerically computed values of the constants are given in Table 1 for up to L = 9. 

The interaction H[ 2 between two closed chains is mediated by L bonds of strength 
J', H' l2 = J' Hi Si,i-Si,2, which in terms of the effective degrees of freedom is obtained by 
J2i=Ai : i. Now it is strightforward to write the effective Hamiltonian for the full system 
within the closed chain ground state subspace as 

Heff = K L J2 S t .S l+1 (b L + + Xfxf + i), (5) 

i=i 

where K L = J'a L L/2, and b L = 2 7 £ (K 3 = 8 J'/3, 63 = 1/8, and K 5 w 2.845J', b 5 w 0.06528, 
see Table 1). It is interesting to note that only xy-component of the chiral interaction enters 
the effective Hamiltonian leaving only a two-fold up-down symmetry for the chiral variables. 
The rotational invariant for the total spin is expected as we started with a fully rotational 
invariant interaction and we did not break the symmetry by restricting the Hilbert space 
to the spin-1/2 sector of every closed chain. Just by looking at the form of the effective 
interaction we learn that on a pair of neighbouring sites the spins and chiral variables are 
strongly correlated. If we try to lower the spin energy by forcing a singlet state between a 
pair of neighbouring sites, the chiral variables go into a triplet state (with x z = 0). However, 
because of the additive constant bi in front of the chirality interaction, spin triplet and a 
chiral singlet is not favored. 

Let us try some variational ansatz. A simple variational state of spin singlets and 
chiral triplets on nearest neighbours, ipsr = Yliodd^i *W wnere V'f is a s P m singlet between 
sites i and i + 1 and ipj is a chiral triplet state with x z — 0, gives an upper bound on the 
ground state energy of the effective Hamiltonian given above E e ff G /K L < E S t{2)/K l = 
-f (6l + |)f (=-0.9375AV4, for L = 3). The factor -3/4 is the spin singlet energy, one half 
is the chiral triplet energy and N/2 is due to the fact that only half the bonds have spin 
singlets and chiral triplets. This state is of course too simple to give us any insight about the 
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structure of the actual ground state. However, as we will see below, it is doing better than 
expected. Let us use Jordon-Wigner fermion language for the chirality terms. The xy-term 
translates into a simple spinless fermion hopping. Now it is tempting to try a free fermion 
ground state for chirality variables, with an average value of - for the xy chiral interaction. 
In this case, we would be left with a Bethe ansatz ground state for the spins, giving an 
energy estimate E BF /K L = -(log 2 - \)(b L + ±)iV(« -0.786iV/4 for L = 3). Compared 
to the singlet-triplet valnece bond variational state the Bethe Free- Fermion variational state 
has 20% larger energy, indicating that a correlated state with spin singlets and chiral triplets 
with longer-ranged valence bonds would do much better. Let us consider a four-site system 
and construct valence-bond state space. We form two orthogonal states for a four-site system 
(the sector with S = has only two states in this case) from the valence bond states, viz. 
A = (12)(34), B = {(23)(41) - i(12)(34)}2/ v / 3, where (12) stands for a spin singlet between 
the two sites. Similarly we construct the chiral states A' and B' where (12) denotes a chiral 
triplet with \ z — 0- Diagonalizing the effective Hamiltonian given in Eq. |5], with this four-fold 
spin-chiral valence bond subspace yields a variational energy of Est{^)/Kl ~ — 0.9695iV/4 
for L = 3. This is a 3% lowering of energy compared to the nearest-neighbour valence bond 
state we considered above. The probability amplitude for a valence bond configuration with 
a bond occupied by a spin singlet and chiral triplet is more than the other configurations, 
implying strong correlations between the spin and chiral degrees of freedom. It is interesting 
to note that this further neighbor valence bond variational state has an energy very close 
to the exact ground state energy by less than three percent (the indication from numerical 
diagonalization of the effective Hamiltonian on a finite site system up to N = 10 is that the 
ground state energy tends to -NK L /4 for L = 3, as we discuss below). 

An effective way of constructing a convergent series of upper bounds on the exact 
ground state energy is by calculating the ground state energy E{n) of a finite open chain 
of n — even sites, giving rise a variational inequality Eq/N < E{n)/n. The best bound 
is obtained from the biggest open chain that can be numerically diagonalized. Since the 
number of states per site is 4, we are able diagonalize only up to ten sites. This gives us a 
bound for L = 3 

E GS /NK 3 < £(10)/10 » -0.9932/4 (6) 

From this calculation it seems that the ground state energy is tending towards -NK^/A from 
above. Similarly one can construct a sequence of lower bounds variationally from /—site 
chain energy through Eq/N > E(n)/(n — 1). However, this sequence is seen to be not as 
convergent as the sequence of upper bounds. A better convergent sequence of lower bounds 
can be constructed by varying the bond strengths for a given n— site chain, which becomes 
very complicated for even n — 6, and we will not present the details here. The ground state 
is seen to be in a sector with S = and x 2 = 0, in all the finite chain diagonalizations we 
have done up to N = 10. The spin gap, defined as the difference between the energies of 
the lowest energy states of the sector S = 1, x z — 0, and S = 0, x z — 0, is seen to fall of as 
the number of sites is increased, as shown in Figure 1, for different values of L. However, 
the data at hand is insufficient to fit the finite- N scaling of the gap. The gap seems to be 
non-monotonic as a function of L for a given value of N, as it can be seen from the figure 
for L = 9 the gap is larger than for L = 3. This could be an artifact of our approximation, 
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implying the perturbation theory in conjunction with the truncation we have implemented 
may not be valid for larger L. We have studied the spin-spin correlation function also for the 
ten-site chain and is seen to be slowly decaying as in the case of L = 1. We may conjecture 
that the decay is likely to be slower for a larger L than that for the case of L = 1, though 
we are not able to substantiate it just now. This is based on that fact that we have a 
long-ranged valence bond ground state, and the probability amplitude for a configuration 
with a valence bond between two sites far apart would be much smaller for a larger L than 
for the case of L = 1. But the characterization of the long-distance behaviour needs an 
investigation of larger chains. The density matrix renormalization approach [?] can be used 
for this; currently an investigation is in progress incorporating this scheme for the finite-size 
spectrum and thermodynamics of larger chains (N > 10) for the effective Hamiltonian. 

Now we consider the case of a non-zero magnetic field in the Hamiltonian given in 
Eq. |]. It is easy to see that the effective Hamiltonian given in Eq. [| will have an additional 
piece, viz. —hJ2 Sf- Let us focus on the case of L = 3. The effective Hamiltonian is valid in 
the regime h < 3 J/4, for L = 3, to ensure that the ground state subspace be composed of 
S z = ±1/2 states of the L— site chains. Now, if the magnetic field is strong enough to polarize 
the spins in the effective Hamiltonian, the ground state would be just a polarized spins and 
a free fermion state for the chiral degrees of freedom. The range of h for which this state has 
lower energy than S = state works out to be for 0.75 J > h > (7n — 8)J'/6n ~ 0.74225 J'. 
The free fermion chiral state has low-lying excited states with a gap of A = (77r — 8)K^/4:N. 
This should be contrasted with the case of L = 1, where there are no low-lying states if the 
magnetic field is of the order of the exchange coupling constants. 

In conclusion, we have established an effective spin-chiral chain Hamiltonian for a 
system of odd-numbered coupled Heisenberg spin chains. The effective interaction has a 
product form, viz. isotropic Heisenberg for spins and only xy interaction for chiral spins. 
The spin gap studied numerically using the effective Hamiltonian for L = 3, 5, 7, 9 falls off 
as the number of sites i increases. Indications from a finite-size study are that long-ranged 
resonating valence bond variational states (with spin singlets and chiral triplets) are suited 
to study the behaviour of the effective Hamiltonian. 

It is a pleasure to thank A. M. Sengupta and M. Barma for discussions. 
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L 


h/J 


b L 


K L /J' 


3 


3/2 


1/8 


8/3 


5 


1.1180 


0.0697 


2.8685 


7 


1.0489 


0.0436 


3.2765 


9 


0.8612 


0.0230 


4.8275 



Table l:The energy gap of a L— site ring, and the wave-function renormalization constants 
that appear in the effective Hamiltonian given in Eq. [| for various values of L. 



Figure Caption 



Figure 1. The spin gap of an L x N system as a function ofl/N for various number of cou- 
pled chains, L = 3,5, 7, 9. The gap is calculated numerically from the effective Hamiltonian 
given in Eq. |5|. 
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